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Abstract 

By integrating the pressure equation at the surface of a self coupled curvi- 
linear boundary, one may obtain asymptotic estimates of energy shifts, which 
is especially useful in lattice QCD studies of nonrelativistic bound states. 
Energy shift expressions are found for periodic (antiperiodic) boundary con- 
ditions on antipodal points, which require Neumann (Dirichlet) boundary 
conditions for even parity states and Dirichlet (Neumann) boundary condi- 
tions for odd parity states. It is found that averaging over periodic and 
antiperiodic boundary conditions is an effective way of removing the asymp- 
totic energy shifts from the boundary. Asymptotic energy shifts from boxes 
with self coupled walls are also considered and shown to be effectively antipo- 
dal. The energy shift equations are illustrated by the solution of the bounded 
harmonic oscillator and hydrogen atoms. 
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I. INTRODUCTION 



In investigating the effects of self coupled boundaries in Monte Carlo lattice QCD simula- 
tions there is a need for guidance as to the expected size and direction of boundary-induced 
energy shifts as well as strategies to remove these effects. Starting with parameterizations of 
(quenched) lattice potentials within the context of a nonrelativistic two-body Schrodinger- 
Pauli Hamiltonian with relativistic corrections, a recent approach uses these potentials to 
produce genuine predictions of charmonium and bottomonium energy levels These 
predictions are obtained from numerical solutions, using a purely spatial grid, of the wave- 
functions plus perturbation theory. Using this approach, one may begin to understand finite 
size effects numerically by changing the size of the enclosing box, which is relatively easy to 
do, unlike the Monte Carlo simulations B. Although the general case must still be studied 
numerically, this article shows that the asymptotic energy shifts from distant boundaries can 
be determined analytically in terms of the unbounded wavef unctions. Where these wave- 
functions are unavailable, it provides functional forms to compare with in potential type 
simulations, and ultimately, directly to the lattice data once box sizes can confidently be 
said to be in the asymptotic region. 

The systems studied here will be described by the Schrodinger equation in the interior 
of a volume, the surface of which is determined by a single parameter of a general sepa- 
rable curvilinear coordinate system. It will be assumed that the fields on the boundaries 
are self-coupled in the sense of being there being either periodic or antiperiodic boundary 
conditions relating the antipodal points at r and — r. We will see that periodic (antiperi- 
odic) boundary conditions require Neumann (Dirichlet) boundary conditions for even parity 
states and Dirichlet (Neumann) boundary conditions for odd parity states for parity invari- 
ant potentials. These boundary conditions are in general not the usual ones applied in cubic 
geometries, which couple opposite walls rather than antipodal points. However, we will see 
that the asymptotic shifts when opposite walls are coupled are effectively antipodal. Thus 
the description here should apply to the situation of Refs. [0,0, which studies confining 
solutions (linear plus Coulomb) in a periodic box. 

We will proceed to a derivation of the energy shift formulas in Sections fTI] and [TTI|. The 
basic energy shift equations for Dirichlet and Neumann boundary conditions will be derived 
in Section [TV]. An asymptotic form of the wavefunction will be presented in Section |V], 
which will then allow us to integrate the differential forms, resulting in simple equations for 
the energies themselves. It will be seen that the asymptotic energy shifts for Dirichlet and 
Neumann cases are equal in magnitude but opposite in direction. In Section we will make 
some simplifying observations regarding surface integrals in the energy shift formulas and 
in Section |V1J| we will explicitly consider the usual lattice situation of cubic geometry with 
self coupled walls. As examples of the use of these formulas, in Section |VIII| we will study 
a one dimensional example, the harmonic oscillator, and a three dimensional example, the 
hydrogen atom. We will close with some general observations regarding qualitative behaviors 
of energy shifts in the general well as some summary comments. 
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II. SCRODINGER MOMENTUM TENSOR DERIVATION 



One way to relate boundary induced energy shifts to properties of the Schrodinger wave- 
functions is to use the momentum tensor. The Lagrangian desity for the Schrodinger case 
may be written 

£ = ^-9^*8^ + (V- E)iP*ip. (2.1) 
2m 

The energy, E, can be viewed as a Lagrange multiplier associated with the normalization 
constraint, / d 3 xip*ip = 1. The Schrodinger equation is recovered from the action defined 
by Eq. ([2.1|) by variation of ip*. (In the time dependent case one concludes ip* and ip are 
related by complex conjugation only after independent variations on ip* and ip.) The mass 
parameter m should be understood to represent the reduced mass in applications to two 
body systems. 

The momentum tensor for complex fields may be written as 

r « = ^^ - + ^* r -^ (2 - 2) 

which is explicitly symmetric in i and j. Using Eq. (|2.1|) , this gives 

Tij = l-ldrtdrf* + d^dnj;*] - ^-d k (rdkij). (2.3) 
As a check, we notice that 

diTij = -iP'djVip, (2.4) 

which is just a form of Ehrenfest's theorem when integrated over the volume. The continuity 
equation 

§f + djTjt = 0, (2.5) 
is satisfied in the time dependent case, where the momentum density is given by 

ffi = ~W-W*), (2-6) 

assuring momentum conservation. We will calculate energy shifts by evaluating the pressure 
exerted by the wavefunction on the surface. For this purpose consider the momentum tensor 
with locally normal indices: 

T nn = r^-0„W - ^~d T ^d T r - t^VV (2.7) 
2m 2m 2m 



Here 



is a locally normal derivative (n is outwardly directed), h n is a possible scale factor, and dr 
are the tranverse derivatives. In the context of integrations over a closed surface, 

_L*.vV + - (2.9) 



giving 



T„„ = ^larfl 2 - (2.10) 
Antipodal boundary conditions are given by 

^(r)| s = ±^(-r)U (2.11) 
^(r)|, = TW(-r)|„ (2.12) 

where the top signs for the periodic case and the bottom signs give the antiperiodic case. 
(One plugs in — r after the normal derivative is taken in ( |2.12j ).) For a parity invariant 
potential, one has 

*l>\. = 0, (2.13) 

by continuity of the wavefunction at the surface for odd parity states, and 

d n ip\ s = 0, (2.14) 

by continuity of the first normal derivative of the wavefunction for even parity states in the 
periodic case. For antiperiodic boundary conditions, the conditions (|2.13| ) and fl2.14|) apply 
instead to even and odd parity states, respectively. When the states can not be classified 
by parity, as when an otherwise good parity system is not centered at the origin, Eqs.( |2.lID 
and ( 2.12 ) can no longer be simplified. 

For the Dirichlet condition Eq. (|2.13|) , one then has 

T nn = ^-\d n ^\\ (2.15) 
2m 

and 

T nn = -T^-rdfy, (2-16) 
2m 

for Neumann conditions, Eq. (|2.14|) . 

Finally, we relate the energy shift to an integration over T nn via 

r r a o r 

AE = P(ao) = J PdV = J da J dsh n P(a), (2.17) 

where 

P(a) = -T nn {a). (2.18) 
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Thus, in a differential sense one has 



da 

resulting in 



-T = - I dsh n T nn (a), (2.19) 



Ta = ~ 1 Li (dsh M\\ (2-20) 



in the Dirichlet case and 



d£ J_ 
da 2m 



dsK^&J), (2.21) 



in the Neumann. Again, in the general case without good parity, the form from Eq.( j2.10 ) 
can not be simplified. 

These differential energy shift formulas are actually exact. The wavefunction, ip, is a 
shorthand for ip(/3(a))\ s , denoting the exact wavefunction on the boundary determined at 
finite /3, which is itself determined by the boundary parameter, a, for a given surface. As 
they stand, these equations are not yet useful because we have no information on the values 
of the wavefunction and it's normal derivatives on S. A hint on how to proceed is provided 
by an alternate derivation of these results, which will now be presented. 



III. PARAMETER VARIATION DERIVATION 

Let us consider the variation of the Schrodinger equation, 

- — V 2 V> + VV> = -Ety, (3.1) 
2m 

with respect to the parameter (not the coordinate), a, which determines the self coupling 
surface, assuming the normalization condition, J d 3 xip*ip = 1, is satisfied. (Usually the 
following procedure determines the normalization integral given E(a), but here we use it in 
reverse.) We have 

- 7T- V> V — + ip V— = —ip ijj + Eip — , (3.2) 
2m oa oa da oa 



after variation followed by ip* multiplication. On the other hand, complex conjugation of 

da 



Eq. (|3.1|) followed by multiplication by |p results in 



J_^VV + £W = E^'. (3.3) 
2m oa oa oa 



The difference produces a perfect differential form, which gives 



~t~ = 7^ / ds I -z-d n ip - V> d n — , 3.4 
da 2m J \oa oa I 
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where d n are again local normal derivatives. (The order of the partial derivatives in the last 
term is immaterial since a appears only in tp.) We then have 

f = l[ d A r , (3.5) 
da 2m J oa 

for Dirichlet boundary conditions, and 

f- = ~ Idsr^rd^, (3.6) 
da 2m J oa 

for the Neumann case. These results may be reconciled with Eqs. (|2.20|) and (|2.21| ) by the use 
of expansions in the normal curvilinear coordinate, n, about the boundary. In the Dirichlet 
case since the wavefunction vanishes on the surface one must have 

^(n-a)|^| s , (3.7) 

resulting in 

^■\ s = -h n d n tfj\ s . (3.8) 

In the Neumann case since the first normal derivative vanishes one has 

d 

d n il) « (n- a)—d n il)\ s , (3.9) 

giving 

^dn^\s = -Kd 2 n ijj\ s . (3.10) 

This reconciles Eqs. (|3.5| ) and ( p.6|) with (|2.20|) and ( |2.21|) . In the general case with no good 
parity, one must proceed differently. One may expand about the boundary the quantity, 

/ d n ip(r) d n ip(-r)\ 



(for \i/j\ s nonvanishing) which also vanishes on S. Then one learns that Eq. (|2.19|) , with the 
general T nn from ( [2.10|) , and Eq.( p.4|) are equivalent when the antipode terms are grouped 
together inside the surface integral. One also learns that one may replace 



(3 - U) 



anywhere in these expressions. 
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IV. BASIC FORMULAS 



The use of the a parameter allows one to derive a differential equation for the energy 
shift, (3, and simplifies the Neumann expression. We begin with a change of variables from 
a to /3(a): 

dip = d(3{a) dip 
da da d/3' 

In the Dirichlet case one may expand the boundary condition, which determines the energy 
values, in f3, 

= ^|.«^(0)|. + /?|||., (4.2) 
(■0(0) specifies the unbounded wavefunction) allowing the approximate replacement 

gi. « (4.3) 

When ( |3.8| ) , ( |4.1| ) and ( |4.3| ) are are used in ( |3.5| ) , the variables j3 and a separate and we have 

PM dp r a o da 

~w K ~ m ] jdsh^imr { ' } 

This identifies the anayltic form for l//3(ao) up to a constant. Thus, asymptotically one has 

/3> ) - ( 2m r rj ^ d * 9 ) " . (4.5) 



\ ■ Jdsh-^(0)\ 2 , 
The Neumann case proceeds similarly, beginning with 



d d/3(a) d 



The boundary condition is again expanded in /3, 



resulting in 



= d n i,\ s w 9 n ^(0)| s + P-Hpd n il>\ 8 , (4.7) 

^d n ^\ s ^~d n m\s. (4.8) 

As we will see in the next Section under very general conditions for good parity states, 

0(/3(a))| s «20(O)| s . (4.9) 
When (|4.6| ), (|4.8|) and (|4.9|) are used in ( |3.6|) , no integration is needed and we simply obtain 

P N (ao) ~ - / ^o0*(O)9 n 0(O), (4.10) 

for the energy shift. Note that we have not assummed that the wavefunctions are separable. 
For a general confining ellipsoidal surface it can be shown that in the asymptotic limit the 
appropriate scale factor goes to one 0. In the following we will also see that it is always 
possible to choose surfaces for which these factors are effectively unity. Therefore we will 
set h n = 1 in the following. 
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V. ALTERNATE APPROACH 



An alternate approach determines the energy shifts from the form of the asypmtotic 
wavefunction. Assuming that the potential, V, is spherically symmetric far from the region 
where the unbounded wavefunction exists (measured by it's average radius, < r >, say), we 
have the approximate wave equation 

W(<0) = 0, (5.1) 
which is solved by 

^{a)) = m\s{e- J > {r ' W ±J> {r ' )dr ') 1 (5.2) 

where r a = r\ s and where the upper sign goes with Neumann boundary conditions and the 
lower goes with Dirichlet and the second term gives a vanishingly small contribution far from 
the boundary. In the immediate vincinity of the boundary, the unbounded wavefunction may 
be characterized as 

m = m\se-^ K(r ' )dr '. (5.3) 
Thus we conclude in the Dirichlet case, 

W03(o))| ( «2W(0)|„ (5.4) 
giving the alternate expression, 

(3 D (a ) « — r da [ ds |<9„^| 2 • (5.5) 
Similarly, in the Neumann case one has (|4.9|) above as well as 

d^(P(a))\ s ^2d 2 n m\s, (5.6) 

giving the alternate form 

2 r a o t 

p N (a ) w — / da dsip*(0)d?Tp(0). (5.7) 
m Joo J 

The next Section will show that use of the aymptotic forms of the wavefunctions allows the 
Dirichlet or Neumann energy shift formulas to be integrated and results in 1 

f3 D (a ) w -/3 N (ao) « ~J ds ^*(0)5 n V(0). (5.8) 

This equation can also be interpreted as the energy shift in systems in higher or lower 
dimensions as long as the basic asymptotic wave equation holds in the appropriate 

variable. 



d 2 
dr z 
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VI. SURFACE INTEGRALS 



In the last two Sections we have derived differently appearing formulas for the energy 
shifts. These will be reconciled in this Section. In the process we will learn about the 
surface integrals in these formulas in the general situation when the enclosing surface is not 
spherical. 

The reconcilation of ( 4.10 ) and (5/7) in the Neumann case is accomplished by use of the 



asymptotic solutions to Eq.( |5.1| ). The condition r a ^> < r > in each quantum state (or the 
appropriate statement in one dimension) and the existence of the asymptotic wave equation, 
Eq.(|5.1|), insures that we can ignore any r dependence in the unbounded wavefunction, ^(0), 
when integrating or taking normal partial derivatives, other than from the exponential in 
Eq. fl5.3p . Writing the general asymptotic form of the unbounded wavefunction (still not 
assuming separability; Q represents the angular variables) as 

tjj(0) = N(r,n)e-lo dr ' K{r '\ (6.1) 
one can show that either (|4.10|) or (|5.7|) give the common form, 



m 



(5 N {a Q ) « / ds K,(r a ) d n r \ip(0)f, (6.2) 



This form may be further reduced because of the presence of the exponential in the surface 
integral. This will suppress contributions from regions of the surface not in the immediate 
vicinity of the points of closest approach to the potential center. Let us call r± this closest 
distance. Thus, we may pull the r dependent quantities outside of the spatial integral, frozen 
at their values for r = r±, 

(3 N (a ) » -- k{t ± ) ?LL / ds \m\ 2 - (6-3) 
m on J 

In the Dirichlet case, Eq. (|5.5| ) may similarly be reduced to the negative of the right hand 
side of Eq. fl6.3|) using the asymptotic form, Eq.fl6.1|). On the other hand, starting from fl4.5|) 
one may write 

ds\m\ 2 = e- 2 ^ K{r ' )dr ' fds\N(r a ,n)\ 2 e- 2l ^ Kir ' )dr '. (6.4) 



The surface integral in ( |6.4|) can be treated as a constant in the boundary integral in ( |4.5| ) 
(verified below), again resulting in the negative of the right hand side of Eq. fl6.3| ) for (3 d {(Iq). 

Eq. (|6.3j ) is still unnecessarily complicated. Since the near points on the boundary will 
dominate the surface integral, there are three cases that can occur for the locus of these 
points: sphere, circle and point. In the sphere case one simply has 

PlhM « -- «(r±) / ds \m | 2 . (6.5) 

In the case where the near points are a circle, one can construct a local cylindrical surface 
for the integration, giving 



9 



J ds\N{r a ,Q)\ 2 e 2 w r ± J J d<p J°° dz \N(r ± , Q±)\ 2 e' 2 ^'^ , 

nr ± #|iV(r ± ,fi ± )| 2 / dze ~ z , 

JO J-oo 

d<p\N(r ± ,Q ± )\ 2 . (6.6) 

Jo 



nr 3 r 2-K 



This results in 

Plde(ao) ~ -- V^i^(^) r d<j>\^ ± (0)\ 2 , (6.7) 



where ?/>_l(0) is the wavefunction along the r± points. Likewise, if the nearest equidistant 
points are discrete, one can consider a locally flat surface and evaluate 

J ds\N(r a ,n)\ 2 e~ 2j rl K{r ' )dr ' ^27r£ |^(r ± , n ± )\ 2 J™ dppe' 2 ^'^, 

^2nJ2\^(r ± ,n ± )\ 2 dppe'-" , 

K>± £ mr ± ,n ± )\ 2 , (6.8) 



where the sum is over the equidistant points. This results in 

/S^o)--^EhMo)l 2 - (6.9) 



This last case is the appropriate one for a cubic box enclosure with antipodal boundary 
conditions. 



VII. BOXES WITH SELF COUPLED SIDES 

The antipodal boundary conditions considered here are not the ones used in Monte Carlo 
lattice simulations since they distinguish one spatial point as the "center" . These simulations 
generally use periodic or antiperiodic couplings on opposite box sides. Nevertheless, the 
antipodal conditions are effectively implemented in such simulations because of the local 
nature of the surface integrals. Let us see how this comes about. 

In general one should consider both continuity of ip and d n ip on the surfaces in construct- 
ing the energy shift formulas for boxes with self coupled sides. Let us consider the energy 
shift from the x = ±a surfaces of the box for a parity invariant, centered potential. The 
asymptotic wavefunction with antipodal boundary conditions in this case can be written 

W(a)) « 1>(r a , +) e -£ K(r ' )dr ' ± V(r a , -) e£ K{r ' )dr ' , (7.1) 
where the upper sign is for periodic and the lower is for the antiperiodic case, and where 

^(r a ,+) = N(r a ,n) e -C^)*', (7.2) 
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gives the general form of the unbounded wavefunction with r replaced by r a = ^/ d 2 + y 2 + z 2 . 
■0( r a> — ) is the same as above except with Q reflected through the x = plane. The second 
term on the right of Eq. (|7.1|) will be exponentially suppressed in the interior of the box. Now 
using the general results, Eqs. (|2.10|) and (|2.19|) , and integrating over a using the asymptotic 
form ( |7. 1| ) , we have 

(3 box (a ) ~ =F— / ds K(r a )d n rRe(r(r a ,+)Hr a ,-)). (7.3) 
m J 

(So indicates both sides of the box.) We will get the full energy shift by integrating over the 
other sides of the box in a similar manner. Because of the local nature of these integrations, 
the contributing portions of the box sides are almost antipodal, and the factors n(r a ) and 
d n r can come outside of the surface integral. Then, for wavefunctions with good parity, 

V±(r a ,+) = ±V±(r a) -), (7.4) 

Eq.([7.3|) in the periodic case just gives the results of the last Section for Neumann and 
Dirichlet boundary conditions, respectively. Antiperidicity just reverses these results. Near 
the box edges the representation ( |7.1| ) fails to be correct, but the local nature of the surface 
integrations allows it as an approximation. When the potential is not centered or the states 
do not have good parity, similar considerations show that periodic and antiperiodic energy 
shifts are negative and equal, at least in the case where one point is much closer to the force 
center than the other. 



VIII. EXAMPLES 

We will consider two examples of the use of the energy shift formulas, the one dimensional 
harmonic oscillator and the hydrogen atom. 

For the harmonic oscillator, the asymptotic form of the wavefunctions are given by 
(x = (mw)" 1 / 2 ; x = Xo/x; Xq is the boundary distance from the origin)) 

~M _I~2 

$ N {x) = (8.1) 

Use of ( |8~7[| ) in Eq. ( |5.8| ) (the "surface" integral simply supplies a factor of two) gives imme- 
diately that 

PSM - [ ) -L, r ■ (8 ' 2) 

On the other hand, the explicit energy shifts for the harmonic oscillator may be found 
from the general forms of the solutions, proportional to confluent hypergeometic functions, 
using the boundary condition on the wavefunction. The form of the wavefunction is Q 

ifj even (x) ~ e~^ 2 F( ~ A 4 +1 \^\x 2 ), (8.3) 
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for the even parity states, and the energies are E = X/2mx 2 . Expanding A as 

A = (4iV + 1) + 2mx 2 AE, 
) one has in the Dirichlet case , ip even (x) = 0, that 



(JV=0,1,2, 



F(-N 


i 

2 


X 


2 ) 




2mx 2 ^F(-N 




A 
4 


1 

2 


5 2 )U=o 



5.4) 



•5) 



Asymptotically, 



F(-iV|-|5 2 ) 



(I) 



(8.6) 



JV 



where 



2V-1 



(a) 



n(a+^) 



T(iV + a) 
Tfa) 



(8.7) 



The method of taking derivatives of confluent hypergeometric functions in Ref. p[ is used, 



neglecting polynomials in xq up to order 2N in ^F(—N — 



A 1 1-2 



4' 2' 2 



I A=0- 



the summand and considerable reduction, this can be put into the form 



Ai#(*o) - f| 



~2r-4Af+l 



v 



After a shift in 



(8.8) 



In dropping the polynomial contribution, we must have that the major contribution to the 

2' 



infinite sum in (|8.8|) occurs for r ^> 2N. Since < N\x 2 \N >= x 2 (N + ^) in the unbounded 



state N and because the major contribution in the infinite sum actually occurs for r rs x 2 , 
we have that ^> 2 < N\x 2 \N > for these asymptotic forms to hold, as would be expected. 

One notices that the Eq. fl8.2p is odd in xq, while QS.8| ) is even, so a term by term com- 
parison of the two expressions (expanding e~ x = (e x ) _1 in a power series) is not possible. 
However, the asymptotic limit of the two sums are identical. This can be established by 
using 



(N + ^) t - 2N (t - 2N) N+1 
x 2t+1 



-i2 



N 



E 



i- 



(2JV-1)!! 



ire 



(8.9) 
8.10) 



r=2JV+l \2' T + 1 



Eq. (|8.9|) holds for r 2N and Eq. (|S.10|) can be justified by using Eq.( |8.7|) , replacing the 
sum on r with an integral, and doing the shift r — > r — | (neglecting a polynomial of order 



AN + 1 in x). Finally, using (8.9) and ( |8.10Q in (8.2), we now have that 

AE%(x )=P° N {x ), 



(8.11) 
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as is appropriate for the even parity states. Going back to the general form of the wavefunc- 
tion, one may also verify that asymptotically AE^(xq) = —AE^(xq) from the Neumann 
boundary condition, d x ip even (x)\ x=xo = 0. 

For our second example we consider a hydrogen atom bounded by a sphere. The asymp- 
totic form of the radial wavefunctions for the hydrogen atom are (f = ro/Nao] 'W is the 
Bohr radius) 



3/2 (of\N-l -f 

Rnl{t) « ( ) i ( 8 - 12 ) 
' '2N{N + L)\{N - L - 1)! 



Na ( 



This immediately implies the Dirichlet energy shifts 

6 D (r ) ~ 1 {2 " )2Ne ' 2f ' (8 13) 

from Eq.( |5.8j ). Obviously, the surface integral was trivial in this case. 

The explicit form of the energy shifts for the hydrogen atom have previously been eval- 
uated in Ref. in the Dirichlet case. This was done by expanding the general form of the 
radial wavefunction 

R(r) ~ e~^r L F(L - n + 1 \2L + 21 — ), (8.14) 

with energy 

E = -(2ma 2 n 2 )- 1 , (8.15) 

to first order in a Taylor series in the energy shift, 

nw JV + ma 2 N 3 AE, (8.16) 

(N = 1,2,3,...) and then solving for AE® L with boundary condition R(r ) = 0. One 
obtains after reductions, 



=27V- 



_ t (r — L — N — l)\(r + L — N)\ 



[8.17) 



using the asymptotic form of the wavefunctions and neglecting a finite polynomial series of 
order iV — L — 1 in r in the denominator. For r ^> N or L (guaranteed for r 3> N 2 ), one 
has 

< T - 2iV - 1 ) ! - 1 (8.18) 



(r-L-JV-l)!(r + L-JV)! rV 
which results in 

A^ L (r )=^ L (r ), (8.19) 

when another polynomial of order 2N in r is neglected in forming the exponential, e 2r . 
Once again, one may verify that the Neumann shift is just the negative of the Dirichlet one. 
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IX. COMMENTS AND SUMMARY 



A number of comments may be made about energy shifts in general from the above 
formulas. 

• Averaging over Dirichlet and Neumann boundary conditions by averaging over periodic 
and antiperiodic boundary conditions in lattice simulations, is an effective way of removing 
the asymptotic energy shifts for nonrelativistic bound systems with good parity. This point 
can and should be tested numerically in potential simulations. It is the assumption of good 
parity which allows this statement to be made, and is not true for antipodal boundary 
conditions in general. 

• The exact formula, Eq. Q2.20| ), shows that the energy shift is never negative in the case 
of Dirichlet boundary conditions, which corresponds to the usual meaning of a bounded 
system. This makes sense from the point of view of the the extra kinetic energy generated 
from the Heisenberg uncertainty principle. 

• The exact differential energy shift in the case of Neumann boundary conditions for a 
central potential bounded by a sphere can be written as 



V{a) _ E {a) + L{L 4 1 ' 



2ma 2 



(9.1) 



Assuming the right hand side is dominated by the (positive) spherical confinement potential, 
V(a), at large a, one immediately sees that the asymptotic energy shifts are negative as a 
is decreased. For purely Coulombic systems the right hand side is instead dominated by 
E(a), which is negative, leading again to asymptotic negative energy shifts. At extremely 
small a, one expects the positive Neumann box energies, E ~ 1/a 2 to dominate, leading to 
positive energy shifts. Thus, in the simplest scenario one expects a minimum in the energy 
at some intermediate value of a. This seems to be what is seen in the numerical solutions 
corresponding to Neumann boundary conditions (periodic even parity states) in Ref. [0, 
although the confining geometry is a cube, not a sphere. Actual proof of this scenario in 
spherical and other geometries however is difficult and more complicated types of behaviors 
can not be ruled out. 

As illustrations of the derived formulas we considered two examples of centered potentials: 
a one dimensional harmonic oscillator and a spherically bounded three dimensional hydrogen 
atom. The asymptotic energy shifts were calculated by the use of the derived formulas and 
verified from the known forms of the wavefuctions for the even parity oscillator and general 
quantum number Coulomb systems. Although we have not illustrated the more general case 
of a confining nonspherical boundary, formulas have been given for the asymptotic surface 
integrals encountered. The common situation of a confining box with self coupled sides is 



also amenable to numerical calculation using a formula of Section |VII| . Thus, the formalism 
here may readily be applied to lattice simulations of nonrelativistic bound systems in the 
usual box geometry. When applied to truly periodic systems, not self coupled as in lattice 
studies, these results actually estimate the location of one edge of the Bloch energy bands 
which emerge ||. 
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FOOTNOTE 

1. An alternate unintegrated expression, coming from the Dirichlet expression, Eq.( |4.5| ), 
is 



a o da 



— 1 
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